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We discuss results of the k, - factorization approach for heavy quark-heavy antiquark correlations 
in proton-proton and proton-antiproton collisions for RHIC, Tevatron and LHC. We consider cor- 
relations in the azimuthal angle as well as in the two-dimensional space of transverse momentum 
of heavy quark and heavy antiquark. We compare results obtained with the help of different 
unintegrated parton distributions (UPDF) from the literature. 
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1. Introduction 

The heavy quark-antiquark hadroproduction is known as very important test of conventional 
gluon distributions within a standard factorization approach. Standard coUinear approach does not 
include transverse momenta of initial gluons. The method to include transverse momenta is kt - 
factorization approach [[l], Q At the leading order of the coUinear approach the heavy quark 
and heavy antiquark are produced back-to-back. In the unintegrated parton distributions (UPDF) 
approach [Q ||, the azimuthal angle and pt decorrelations (from the coUinear leading-order 
configurations) are obtained already in the leading order of perturbative expansion [0]. In Ref.[^ 
we have explored in detail cc kinematical correlations at the Tevatron energies. Here we wish to 
discuss a more general case of (heavy quark)-(heavy antiquark) correlations for the RHIC, Tevatron 
and LHC energy range. We wish to empasize that this is not yet well explored field of high-energy 
physics which could and should be studied in the future, in particular at LHC. 

2. (Heavy quark)-(heavy antiquark) correlatons 




Figure 1: Typical diagrams for heavy quark - antiquark production. 



Let us consider the reaction hi+h2-^Q + Q+X (see Fig.^, where Q and Q are heavy quark 
and heavy antiquark, respectively. In the leading-order kf- factorization approach the cross section 
in rapidity of Q (ji), in rapidity of Q (j2) and transverse momentum of heavy quark (pi,f) and 
heavy antiquark {p2,t) can be written as 

d(J y r d^Kijd^K2^t 1 . 1 2 

dy\dy2(P' Pi,td^ P2,t J 7i K l6K'^{x[X2sy '■' 

5^{Kl^t + 1<2,t-Pl,t-P2,t) fi{xuKlj) fj{x2,Klt) , (2.1) 

where fi{xi,Kfj) and /y(.X2,K"|,) are so-called unintegrated parton distributions. The two extra 
factors l/n attached to the integration over d^Ki t and d^K2^t are due to the conventional relation 
between unintegrated and integrated parton distributions. The two-dimensional delta function as- 
sures momentum conservation. 
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The unintegrated parton distributions must be evaluated at: x\ = ^ (exp(3;i) +exp(j2))> X2 = 

V ^ 

^ (exp(— ji) +exp(— j2))- The matrix element must be calculated for initial off-shell partons. 
The corresponding formulae for initial gluons were calculated in [|l], ^ (see also [|3|]). In our paper 
we compare results obtained for both on-shell and off-shell matrix elements for charm-anticharm 
correlations. 

3. Results 




Figure 2: da/dcj) for charm - anticharm production at W = 1960 GeV for different UGDFs. The results 
with on-shell kinematics are shown in the left panel and results with off-shell kinematics in the right panel. 
In this calculation both factorization and renormalization scales were fixed for 4 m^. 



In Fig.2 we compare results for different unintegrated gluon distributions from the literature. 
Different results are obtained for different UGDF. BFKL dynamics leads to strong decorrelations 
in azimuthal angle between charm and anticharm quarks. In contrast, the nonperturbative GBW 
glue leads to strong azimuthal correlations between c and c. The saturation idea inspired KL 
distribution leads to an local enhancement for 0cc ~ which is probably due to simplifications 
made in parametrizing the KL UGDF. In the last case there is sizeable difference between the 
result obtained with on-shell (left panel) and off-shell (right panel) matrix elements. All this is due 
to an interplay of the matrix element and the unintegrated gluon distributions. 

In Fig. 3 we show results for correlations in (pij) of c and {p2,t) of c for different uninte- 
grated gluon distributions. Our results depend strongly on different UGDF. This is very interesting, 
because it can be verified in future experimental studies. 

In Fig.4 we show another type of two-dimensional distribution. Generally, the bigger trans- 
verse momentum of the produced quark (or antiquark) the stronger back-to-back correlation in 
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azimuth. Similarly as in the previous case, the results depend on UGDFs used in the calculation. 
The somewhat stronger back-to-back correlation in the case of the Kwiecihski UGDF requires a 
separate discussion. As an example in Fig.5 we show the dependence of the azimuthal correlations 
for pi^t,P2,t £ (5,15) GeV. We show the result for different values of the parameter bo describing 
initial (for QCD evolution) nonperturbative distribution of gluons in nucleons as well as for differ- 
ent values of the scale parameter n^. The results can be characterized as follows. The smaller bo, 
the bigger the decorrelation between c and c. On the other hand, the decorrelation increases with 
raising Physically one would expect = fjL^{pi^t,P2,t)- A reasonable choice would be 

^lHpi^,P2,t)=Ciml + C2{plt+plt) , (3.1) 

where Ci , C2 ~ 1 may be expected. This would require a use of running scale instead a fixed one 
as at present. In the moment this was not possible technically and therefore only typical values of 
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Figure 4: Two dimensional distributions in pt and (j) for KMR, Kwieciriski (bo — 1 GeV ^, fJ.^ — 10 GeV^), 
BFKL, KL for W = 1960 GeV. 



the scales were taken (see Fig.§). The scale fx^ = 100 GeV^ is justified by the averaged values of 
the momenta in our calculation {p\t)^ (Pit) ~ 100 GeV^. The figure shows that QCD evolution 
embedded in the Kwieciriski equations leads to a strong extra decorrelation of c and c in addition 
to that of the nonperturbative origin as encoded in the parameter bo. 

It was predicted for jet-jet correlations some time ago that a large rapidity gap between jets 
leads to increased decorrelation in azimuth. Can one expect a similar interesting effect for cc 
correlations? In Fig.^ we show our result for cc correlations as a function of the rapidity gap 
between c and c. Of course, the bigger rapidity distance between c and c the smaller the cross 
section. However the shape of the distribution in azimuth stays practically the same, except at 
Ay « 0, where a singularity for collinear quark and antiquark emissions shows up. 

Finally we wish to come to azimuthal correlations of t and t. Such an analysis seems possi- 
ble even at present at the Tevatron. In Fig.^ we show the azimuthal correlation function for two 
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Figure 5: (Color on-line) The azimuthal correlations for the gg cc for RHIC and Tevatron energy with 
the Kwiecinski UGDFs for different values of the nonperturbative parameter bo (0,5 GeV^', 1 GeV^', 
2 GeV^i) and for different evolution scales = iq GeV^ (blue line) and = iqo GeV^ (red Hne). The 
initial distributions ( without evolution) are shown for reference by black lines. 



energies: W = 1960 GeV (Tevatron) and W = 14 000 GeV (LHC). We show separately two contri- 
butions: gg tt (red) and qq — > tt (black). While at the Tevatron energy the qq — > tt contribution 
dominates over the gg tt contribution, the situation reverses at the LHC energy where the gg tt 
is the dominant contribution. The shape of both components is, however, very similar. 

It would be very interesting to compare the result obtained with in the /:rfactorization approach 
here with the result of the NLO collinear-factorization approach. The work in this direction is in 
progress. 



4. Summary and outlook 

Summarizing, we have discussed the (heavy quark) - (heavy antiquark) correlations in proton- 
proton and proton-antiproton collisions within the -factorization approach. Different unintegrated 
gluon distributions have been used in the calculation. The results for azimuthal angle as well as in 
P\,t>^ P2,t correlations have been presented. 

In principle, corresponding experimental results would provide new information and could 
test the unintegrated gluon distributions from the literature. In practice, this could be more difficult 
as one measures rather heavy mesons or charged leptons from the decays and only correlations of 
such objects can be studied. However, both heavy mesons and charged leptons from the decays 
"remember" to large extend the direction of initial heavy quark/antiquark. 
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Figure 6: da/d(j) versus rapidity gap between c and c for W = 1960 GeV^. (pi j = p2.t (5, 15) GeV, 
Ay = yi -y2, ji^ = 10 GeV^ (left panel), ji^ = 100 GeV^ (right panel). 
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Figure 7: (Color on-line) da/d^ for tt production at W = 1960 GeV (left panel) and W = 14000 GeV 
(right panel) for two process: gg tt and qq — > tt for Kwiecinski distributions. In this calculation /i^ = 
10000 GeV^. 



On the theoretical side, a relation between /cf -factorization and standard NLO approaches 
should be better understood and clarified. 

We expect exciting future studies of QQ correlations at the LHC, especially for tt where de- 
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tailed studies will be possible for the first time in the history. 

Acknowledgments This work was partially supported by the grant of the Pohsh Ministry of 
Scientific Research and Information Technology number 1 P03B 028 28. 

References 

[1] S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. 366 (1991) 135. 

[2] J.C. Collins and R.K. Ellis, Nucl. Phys. B360 (1991) 3. 

[3] R.D. Ball and R.K. ElUs, J.H.E.R 0105 (2001) 053. 

[4] D. Kharzeev and E. Levin, Phys. Lett. B523 (2001) 79. 

[5] K. Golec-Biemat and M. Wiisthoff, Phys. Rev. D60 (1999) 114023-1. 

[6] M.A. Kimber, A.D. Martin and M.G. Ryskin, Eur. Phys. J. C12 (2000) 655; 
M.A. Kimber, A.D. Martin and M.G. Ryskin, Phys. Rev. D63 (2001) 114027-1. 

[7] M. Luszczak and A. Szczurek, Phys. Rev. D73 (2006) 054028. 



8 



